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This is an extended and pedagogically oriented version of our recent work in ref. [1] 
where we proposed an improvement of the splitting functions at small x which 
overcomes the apparent problems encountered by the BFKL approach. 



1 Introduction 

The theory of scaling violations in deep inelastic scattering is one of the 
most solid consequences of asymptotic freedom and provides a set of fundamen- 
tal tests of QCD. At large imd not too small but fixed x the QCD evolution 
equations for parton densitiesH provide the basic framework for the description 
of scaling violations. The complete splitting functions have been computed in 
perturbation theory at order Ug (LO approximation) and (NLO)lj. For the 
first few moments the anomalous dimensions at order are also knowncl. 

At sufficiently small x the approximation of the splitting functions based 
on the first few terms of the expansion in powers of Ug is not in general a 
good approximation. If not for other reasons, as soon as x is small enough 
that QfsC 1: "^itti C = logl/a;, all terms of order as(as^)" and a^(as^)" 
which are present u in the splitting functions must be considered in order to 
achieve an accuracy up to order a^. In terms of the anomalous dimension 
7(iV, Q!s), defined as the iV-th Mellin moment of the singlet splitting function 
(actually the eigenvector with largest eigenvalue), these terms correspond to 
sequencesjpf the form (ag/N)" or as{as/N)^. In most of the kinematic region 
of HERAQ the condition ~ 1 is indeed true. Considering that, by kine- 
matics, as£, < aslogs/Q^, and s « 10^ GeV^, as{m\) « 0.119, we see that 
at Q2 = 3^ 10^ 102, 103 GeV"^ a,^ can be as large as a^^ < 4.3, 3.0, 1.2, 0.6, 
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respectively. Hence, in principle one could expect tQ,|aee in the data indica- 
tions of important corrections to the approximation of splitting functions 
computed only up to order and the corresponding small x behaviour. In 
reality this appears not to be the case: the data^an be fitted quite well by the 
evolution equations in the NLO approximationBu. An idea of the quality of the 
fit can be obtained from fig. 1 where a comparison of the data with the QCD 
NLO scaling violations is displayed. Of course it may be that some corrections 
exist but they are hidden in a redefinition of the gluon, which is the dominant 
parton density at small x. While the data dQ_not support the presence of 
large corrections in the HERA kinematic region t3 the evaluation of the higher 
order fiorrfictions at small x to the singlet splitting function from the-BEKL 
theory 0ll30 appears to fail. The results of the recent calculationllilOOj of 
the NLO term xi of the BFKL function x — cksXo + ctsXi--- show that the 
expansion is very badly behaved, with the non leading term completely over- 
throwing the main features of the leading term. Taken at face value, these 
results appear to hint at very large corrections to the singlet splitting function 
at small x in the region explored by HERAEZl. 

Here we review this problem and propose a procedure to construct a mean- 
ingful improvement of the singlet splitting function at small x, using the infor- 
mation from the BFKL function. After recalling some basic notions we start 
by defining an alternative expansion for the BFKL function x(Af) which, un- 
like the usual expansion, is well behaved and stable when going from LO to 
NLO, at least for values far from M = 1. This is obtained by adding suitable 
sequences of terms of the form (ag/M)" or as{as/M)"- to asXo or a'^xi re- 
spectively. The coefficients are determined by the known form of the singlet 
anomalous dimension at one and two loops. This amounts to a resummationlia 
of (as log Q^/^^)" terms in the inverse M-Mellin transform space. This way 
of improving x is completely analogous to the usual way of improving 7c3. 
One important point, which is naturally reproduced with good accuracy by 
the above procedure, is the observation that the value of x(M) at M = is 
fixed by momentum conservation to be x(0) — 1. This observation plays a 
crucial role in formulating the novel expansion and explains why the normal 

BFKL expansion is so unstable near M = 0, with xo ~ xi ^/M^ and 

so on. This rather model-independent step is already sufficient to show that 
no catastrophic deviations from the NLO approximation of the evolution equa- 
tions are to be expected. The next step is to use this novel expansion of x to 
determine small x resummation corrections to add to the LO and NLO anoma- 
lous dimensions 7. Defping A as the minimum value of x, x(M,„i„) = A, and 
using the results of ref. Eil, a meaningful expansion for the improved anomalous 
dimension is written down in terms of xoj Xii cind ^- The large negative cor- 
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Figure 1: Next-to-leading order QCD fit to HERA and pre-HERA data on the proton 
structure function F2(x,Q'^) [9]. 
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rection to Aq/os = xo(l/2) induced by asXii that is formally of order as but 
actually is of order one for the relevant values of as, suggests that A should be 
reinterpreted as a nonperturbative parameter. We conclude by showing that 
the very good agreement of the data with the NLO evolution equation can be 
obtained by choosing a small value of A, compatible with zero. 

2 Evolution at Small x 

The behaviour of structure functions at small x is dominated by the singlet 
parton component. Thus we consider the singlet parton density 

G(C, t) = x[g{x, Q2) + kg(g> q{x, Q2)], (1) 

where ^ = log 1/x, t — log Q^Z/i^, g{x, Q^) and q{x, Q^) are the gluon and sin- 
glet quark parton densities, respectively, and kq is such that, for each moment 

/•OO 

G{N,t)= / dee-^«G(e,t), (2) 



the associated anomalous dimension ^{N,as{t)) corresponds to the largest 
eigenvalue in the singlet sector. At large t and fixed ^ the evolution equation 
in A^-moment space is then 

j^G{N,t)^j{N,as{t)) GiN,t), (3) 

where as{t) is the running coupling. The anomalous dimension is completely 
known at one and two loop level: 

-f{N,as) = as"fo{N) + alli{N) + ... . (4) 

As 7(A^, as) is, for each A^, the largest eigenvalue in the singlet sector, momen- 
tum conservation order by order in as implies that 

7(l,a.) = 7o(l)=7i(l) = - = 0. (5) 
The solution of eq.(||) is given by 

r^sit) ^(N a') 
G{N, t) = G{N, 0) exp / da', (6) 

where a, = agiO) and the QCD beta function /3(as) is defined by 

-^ = /3{asm (3{as)^-bal{l + b'as + ; 6=^ = ^ 

at An IzTT 

(7) 



At one loop accuracy 'y{N,as) — as'yo{N) and the solution reduces to: 



G{N, t) = G(iV, 0) e(^°"/^) = G{N, 0) [^]^°/^ 77 = log ^ (8) 

Of formal interest is the solution in the limit of fixed coupling. In general, 
in this limit, we have G{N,t) — G{N,0) exp['y{N, as)t] and, at one loop, 
G{N,t) ~ G{N,0) exp(7oQ!st). Comparing with the running coupling case 
at one loop, we see that agt at fixed coupling is replaced by rj/b at running 
coupling. In fact by expanding we have 

77 = log = log (1 + bad) = ba^t + (9) 

It is useful to keep in mind this translation rule for one loop results. 

We now consider the small x behaviour of the one loop solution. Small 
X means large ^, hence small N. At small N, 7o(iV) — ^ + . . .. The 1/A^ 
behaviour correspond to the 1/x behaviour of the one loop singlet splitting 
function Poix), the Mellin transform of 70: 7(A^) = f^dx x^P{x). In general 
P{x) ^ l/a;(log(l/x))" corresponds to 7(A^) ^ nl/N^^-^. We can add a con- 
stant term to 70 (N) in such a way that the momentum conservation condition 
at = 1 is respected: 

7o(iV)«-(^-l) + (10) 

This addition provides a good approximation in a wide region extending from 
A up to A^ ~ 1. From eq.(H) one obtains the leading behaviour: 



G{N,t) = G{N,0) exp 



TT b^N ' 



(11) 



By taking the inverse Mellin transform we obtain: 



G{^,t)= " p(A^e+^(T7-l))G(^,0) (12) 



C — ZOO 



27rj ' 7r6 ' A 



where the integral is along the imaginary axis at RealN = c. The leading 
behaviour at small x is obtained by the saddle point method (see Appendix). 
Let A'o be the point where the derivative of the exponent vanishes: 
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At A'o the value of the exponent is: 



For G{N, 0) regular at > 0, the resulting small x behaviour is given by: 

G(^, t) = GiNo,0)vk-'^ exp ^^exp (-^) (15) 

This is the well known "double scaling" behaviour, first predicted in ref0 and 
developed in refsu. The exponential factor 

T logQVA^ .-..x 

produces a peak at small x which increases with Q^. The x behaviour is 
weaker than any power of x but stronger than any power of a logarithm. 
This one loop prediction is not much modified by NLO corrections because 
the computed two loop anomalous dimension 71 (iV) also behaves like 1/N. 
And in fact the complete NLO QCD evolution reproduces the double scaling 
exponential rise quite accurately. As already mentioned the HERA data are 
quite well fitted for > 1 — 4 GeV^ by QCD evolution computed at 
NLO accuracy. This agreement of the data with the usual evolution equations 
is however rather unexpected. In fact the evolution takes into account all 
terms of order [a^ log 1 /x log Q^]" and [as log Q^]" and are valid for sufficiently 
large and x small but fixed. However possible terms of order [aslogl/a;]" 
are not included in the splitting function, so that the results are not necessarily 
reliable for sufficiently small x at fixed Q^. The powerful approach based on 
the BFKL equation provides a suitable tool towards an estimate of these terms 
and will now be discussed. 



3 The BFKL Function 

The starting point of the BFKL theory is that, at large ^ and fixed (in 
particular at fixed as), the following x evolution equation for M moments was 
proven valid: 

^G{i,M)=x{M,as) G{^,M), (17) 

where 

/oo 
dte-'"G{tt), (18) 
-00 
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For convergence, we consider < Af < 1 as the physical region and assume that 
G{£,,t) vanishes as in the photoproduction Hmit ^ and approaches 
a constant modulo logs in the limit of large Q^. In eq.(|l7|) x(Af, as) is the 
BFKL function which is now known at NLO accuracy: 

x{M,as) ^ asXo{M) + a^,Xi{M) + . . . . (19) 

The BFKL function has been computed in perturbation theory at one and 
two loops starting from the large s behaviour of the amplitude for gluon-gluon 
scattering. 

We consider T, the absorptive part of the A+B — > A+B forward scattering 
amplitude from gluon exchange, as depicted in fig. 2. A and B are generic real 
or virtual particles. In the case of interest, A will be the virtual photon and 
B the proton. From the optical theorem T is proportional to the total A + B 
cross section a. We are interested in the gluon exchange component which 
is dominant in the singlet sector. So in the deep inelastic scattering case the 
virtual photon A is attached to a quark loop which is included in (j)A- Going 
down vertically then there is the g-g amplitude F^v and finally the proton blob 
(^B. The cross section a can be written as: 

_^,(,,) j I- (^)» r„(..,fe, (20) 

Here ki^2 are transverse momentum vectors of the s-channel gluons with vir- 
tuality — (fc)^, 4>A,B are suitable hadronic structure functions, F^r is the g-g 
kernel and s = 2pAPB is the squared invariant mass of the colliding particles. 
Note that a symmetric scale choice was adopted in the s factor. The g-g kernel 
in the MS scheme obeys the generalised BFKL equation: 

NTNih,k2)^S'^^\ki-k2)+ [ £k K{k^,k)TN{kM) (21) 



The BFKL equation can be cast in the form of the x evolution equa- 
tion (17) inverting the iV-Mellin transform eq. ^ and taking an Af-Mellin 



transform eq. (|18|). The inversion of the A^-Mellin turns multiplication by N 
into differentiation with respect to ln(l/x), and the Af-Mellin turns the con- 
volution with respect to k"^ into an ordinary product. It then follows that the 
x(A'/, as) function is related to K by Af-Mellin transformation: 

x(Af, a.) = j d?k K{h,k) i^f^-^ (22) 
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Figure 2: The process A + B 



A + B. 



At NLO and beyond, x(M, as) can acquire a further dependence on kf due to 
the running of the coupling. From the symmetry under fci k2 exchange of 
r^v, the symmetry of K also follows and in turn, at leading order, this implies 
that Xoi^I) = Xo(l ^ M) We have slightly cheated in this derivation in that 
the variable N in eq. ( |2l| ) is defined by the moment integration of eq. (po|), 
i.e. with the symmetric factor {-^f^)'^ ■ The N moments eq. (||) instead are 
taken with the factor — {-^)~^ where is the virtuality of one of the 
two particles, say = k\. This difference-is irrelevant at LO, but at NLO 
it must be taken into account and modifies llj the relation between x ^^icl K 
eq. (H). 

At one loop accuracy xo{M) is scheme and scale independent and has the 
following simple form (see fig. 4 below): 

Xo(Af) = ^ r + z-^'^ - 2] (23) 

(an alternative expression is Xo(M) = rtc/7r[2-!/'(l) — -0(M) — — M)] with 
tfj{M) = r'(M)/r(M), where F is the Euler gamma function). By expanding 



8 



the denominator and performing the integration, one finds the expression: 



1 1 9 1 

= tI:1mT^+;7TT^— 1 - liJ+2C(3)M=+2C(5)M^+ 

n— 

(24) 

Also note that 

Xo(i) = -41n2 (25) 

TT 



In eq. (17) the coupling as is fixed. The inclusion of running effects in the 
BFKL theory is a delicate point. If we are only interested to next-to- leading 
log 1/x accuracy, the running of the coupling in the evolution equation (^ 
can be expanded out to leading order in t: as{t) — Q;s(/i^)(l — bas{^i^)t). 
Upon M-mellin transformation, this corresponds to a differential operator: 
as{M) = as{l + bag-j^), which can be used in the x evolution equation (p^. 
This affects the relation between the t and x evolution equations, and in par- 
ticular it modifies the relations between x ^^id the anomalous dimension 7 
which we will derive in the next section (duality relation). However, for all 
practical purposes, this modification can be simply viewed as the result of 
an additional contribution to x- once this contribution is incorporated in Xj 
the anomalous dimension 7 is given by the duality relation given in the next 
session. Henceforth, we assume that x includes such contribution, at the ap- 
propriate order in the small x expansion. To next-to-leading order in as (i.e. 
to NLLx), this corresponds to a contribution to xi proportional to the first 
coefficient /?o = — of the /3-function. Since the extra term depends 
on the definition of the gluon density, it is also necessary to specify the choice 
of factorization scheme: here we choose the MS scheme, so that the xi that 
we will consider in the sequel is given by til 

XiiM) = j^nl~5{M) + ^/3one((2^'(l) - ^'(M) - ^'(1 - M)) + j^Xo{M)', 

n ^^^^ 

where the function S is defined in the first of ref. IIJ. Note that beyond the 
leading order the symmetry of the BFKL function under AI ^ 1 — AI is 
destroyed by the effects of running and by the fact that the scales ki and /c2 
are very different in the actual case of interest. As a result the actual BFKL 
function which is relevant for deep inelastic leptoproduction is not perturbative 
near M = 1, the point where the soft scale is approached. 

To start exploring the physical significance of the x evolution equation, we 
work in the fixed coupling approximation. The solution of eq. iRTi) is given by: 



G(e, M) = G(0, M)exp [x(M, a.)C] (27) 
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By taking the inverse Mellin transform, we obtain: 



Jc-ioo 27ri 



xiM,a,)^]G{0,M) (28) 



For t and ^ large we consider the saddle point condition (for simplicity, we 
restrict our discussion to the LO): 

,T7(Mt + asXo{M)A =0, or 1 = -a,Xo(Mo) (29) 



With t and ^ both large we can study three different interesting limits: 1) 
i/^ large and positive (i.e. the limit ^ oo, ^ large but fixed); 2) t/£_ 
(or the Regge limit: a; ^ and large but fixed); 3) t/^ — oo (or the 
photoproduction limit: ^ and ^ large but fixed). In case 1) at Mq the 
derivative must be large and negative, thus Mq is small (see fig. 4). One finds: 



^^0 = Y G(e, t) cc exp [Afot + a.Xo(Afo)C] « exp ^ (30) 

Since we know that agt at fixed corresponds to rj/b for running as, we see 
that the double scahng result, eq.(p^, is reproduced. In case 2), t/(, and 
the derivative must vanish at Mq, or Mq = 1/2. The asymptotic behaviour is 
fixed by xo(l/2) (see eq.(p5|)) and we have: 

G{tt)(xexp[asXail/2)^]~x~^°, Ao=a.— 41n2 (31) 

TT 

This is the hard Pomeron prediction (as opposed to the soft Pomeron, or 
Ao ~ 0), which is well known but of dubious physical relevance, as we shall see 
in the following. In case 3) t/^ is large and negative, so the derivative must be 
large and positive, which leads to Mq « 1 where Xo{M) ~ nc/[7r(l — M)]. A 
simple calculation leads to: 



*^o = l-y^^, ^(CO^Q'expY y (32) 

This result is qualitatively encouraging, because it leads, to the correct linear 
vanishing of G(^,t) in in the photoproduction limit Ej. But, quantitatively 
we cannot trust this result, because non perturbative effects must be important 
at K, 0. Actually this confirms that the BFKL function must contain non 
perturbative effects near M — \ otherwise photoproduction would also be 
perturbative. 
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4 The Duality Relations 

In the region where and 1/x are both large the t and evolution equations, 
i.e. eqs.flllTl), are simultaneously valid_and their mutual consistency requires 
the validity of the "duality" relation E3o: 

X(7(^,«.),a.) = iV, (33) 

and its inverse 

7(x(M,a,),a,) = M. (34) 

In order to derive the duality relations we start from the x evolution equation, 
eq. (|l^) , and take the x Mellin transform of both sides: 

G(e, M) = x(M, «.) G(iV, M) (35) 

Integration by parts leads to: 

[e-^« G(e, M)]o~ + NG(N, M) = «.) G(7V, M) (36) 

The square bracket on the l.h.s is determined by ^ ^ and provides an N 
independent boundary condition, which we denote by —Hq(M) (recall that 
the X evolution equation is only valid at large ^). Solving for G{N,M), we 
find: 

The position of the singularity in N fixes the large t behaviour of G{N, t): 

G{N,t)= / expMtG{N,M) - x{Mp,as) = N (38) 

Jc-ioc 27^^ t-=o 

while the contribution of additional singularities further on the left of Mp is 
suppressed by powers of Q^. Then, at fixed Ug, this must coincide with 

G{N,t) - e'^(^'"^)* (39) 

where 7 is the anomalous dimension function at all order in a^. Thus Mp = 
7(A^, as) and the duality relation eq.(^), and its inverse, eq. (^^ , are obtained. 
These relations are true for the complete leading twist contribution in the do- 
main where the t and ^ evolution equations, i.e. eqs.(|^,[l7|), are simultaneously 
valid. 
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Figure 3: Graphical representation of different expansions of (a) 7 and (b) x ™ powers of 
as and l/N (a) and 1/M (b) to order m and n respectively, and the different relations 
between these expansions. Vertical lines correspond to terms of the same fixed order in 
as: for example the one loop anomalous dimension 70 contains terms with m = 1, n = 
1, 0, —1, —2, .... Diagonal lines correspond to terms of the same order in Os at fixed ^ (a) 
or ^ (b): for example ■ys{as/N) contains terms with m = n = 1, 2, 3, . . .. The sum of terms 
in a vertical line of the 7 plot is related by duality eqs. (^sj^) to the sum of terms in a 
diagonal line in the x plot and conversely (marked by the same line style). The solid lines 
denote terms of the same order in the "envelope" or "double leading" expansion discussed 
in the text. 
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Using eq. (p3[), knowledge of the expansion eq. ( p^ ) of x(M, a^) to LO and 
NLO in as at fixed M determines the coefficients of the expansion of 7(iV, as) 
in powers of as at fixed j^: 

^{N,as)^^s{^)+asjss{^) + ..., (40) 

where 7^ and jss contain respectively sums of all the leading and subleading 
singularities of 7 (see fig. 3). To derive the expressions of 7s and 7ss we start 
from eq. (p3|) written in the form: 

asXobsm + asjssi^) + ■■■] + alxibsi^) + ■■■] + ■■■■ = N (41) 
We then expand in as at as/N fixed: 

N 

Xo(7s(t)) + «sXo(7.(t)) 7..(t) + «.Xi(7.(t)) + = - (42) 



Ois 



From this we find the relations: 

N 



Xo(7s(t)) = (43) 



as 



„, (as\ Xl(7s( N )) (AA\ 



This corresponds to an expansion of the splitting function in logarithms of x: 
if for example we write 



■'.i'i^-Y.A-'iW (45) 

k=l 



then the associated splitting function is given by 

PsiasO^ / ^e^«7.(t)-E7rNv("^^)^'"'^ (46) 

and similarly for the subleading singularities Pss{as£,): etc. From eq.(p4[) it 
follows that 

5^^^ = njvr, g(^) = g^'^ = , g['^ = 2C(3)n,/^, . . . (47) 

We see that the logaritmic terms in Ps only start at order due to the 
vanishing of and g^' . It is also important to observe that the coefficients 
of the Ps expansion are factorially suppressed with respect to those of 7s . This 
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implies that while the convergence radius of the 7s expansion is finite, that of 
Pg is infinite. 

At this point it is interesting to make a digression and compare the small 
X behaviour of the singlet structure functions in the spacelike region with 
that of singlet fragmentation functions in the timelike region. In the leading 
log approximation the timelike splitting functions, relevant for the evolution 
of fragmentation functions, are the same as the spacelike splitting functions. 
Therefore the solution to the one loop evolution equation in the small x region is 
valid both for parton distributions in the spacelike region and for fragmentation 
functions in the timelike region. But if this were the correct asymptotic result 
at x — !■ for the singlet framentation functions it would imply that the layerage 
multiplicity (the N = moment) is singular. Indeed it is in fact knownllB that, 
in the timelike region, the behaviour of the anomalous dimension near N = 
is actually modified by higher order terms according to: 



In X space this is equivalent to the occurrence of terms of order as [as log^ 1/x] "~ ^ 
in nth order perturbation theory in as ■ Thus there are two powers of log for 
each as in the timelike region instead than one. These terms modify the be- 
haviour of j{N,as) near = into a non singular one. This corresponds to 
the well known behaviour of the average multiplicity in gluon jets given by: 



r=(*) 7(0, «.) ^ 2n,loggVA^ 

exp / das —T7 — ^ ~ exp W — ^ITj- (49) 

7a, P{as) y 7r61og^VA2 

which is well supported by experiment. How is it possible that there is such a 
difference between the spacelike and timelike regions? That the fragmentation 
functions are more singular than the structure functions is due to the fact 
that in leptoproduction the final state is totally inclusive. Instead, when we 
compute the fragmentation of a gluon into a gluon, we fix the momentum of 
the observed soft gluon and so its associated infrared singularity cannot be 
canceled against the corresponding virtual contribution and appears in the 
result as an extra power of f . 

Going back to the duality relations, the inverse duality eq. (^4|) relates the 
fixed order expansion eq. (^) of 7(A'^, as) to an expansion of x(-^'^: Qfs) in powers 
of as with ^ fixed: if 

X(M, as) = x.(i7) + asXssm + ■■■, (50) 
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where now Xsijf) and Xssijj) contain the leading and subleading singularities 
respectively of x(M, as), then 

loixsm) = ^, (51) 

We now discuss the physical implications of these formal results for the singlet 
splitting function at small x. 

5 Improving the BFKL Expansion 

In principle, since xo a-nd xi are known, they can be used to construct an 
improvement of the splitting function which includes a summation of leading 
and suMeajding logarithms of x. However, as is now well known, the calcu- 
lation llJ'LaO of xi has shown that this procedure is confronted with serious 
problems. The fixed order expansion eq. ( p^ ) is very badly behaved: at relevant 
values of as the NLO term completely overwhelms the LO term. In particular, 
near M = 0, the behaviour is unstable, with xo ~ ^/M, xi ~ — 1/M^. Also, 
the value of x near the minimum is subject to a large negative NLO correction, 
which turns the minimum into a maximum and can even reverse the sign of x 
at the minimum: in MS for = 3 or 4, we approximately have 

X(i, as) ~ 2.65as(l - 6.2as + ...) (53) 

Finally, if one considers the resulting 7s and 7ss or their Mellin transforms 
Ps{x) and Pss{x) one finds that the NLO terms become much larger thaji the 
LO terms and negative in the region of relevance for the HERA datatZl. We 
now discuss our proposals to deal with all these problems. 

Our first observation is that a much more stable expansion for x(M) can 
be obtained if we make appropriate use of the additional information which is 
contained in the one and two loop anomalous dimensions 70 and 71. Instead 
of trying to improve the fixed order expansion eq. (^) of 7 by all order summa- 
tion of singularities deduced from the fixed order expansion eq. (^) of x, we 
attempt the converse: we improve Xoi^) by adding to it the all order summa- 
tion of singularities Xs CQ- ( |5l| ) deduced from 70, xi(-^) by adding to it Xss 
deduced from 71 eq. (|5^), and so on. It can then be seen that the instability 
at M = of the usual fixed order expansion of x was inevitable: momentum 
conservation for the anomalous dimension, eq. (|^), implies, given the duality 



relation, that the value of xi^^) at Af = is fixed to unity, since from eq. (33) 
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we see that at = 1 



X(7(l, "s), as) = x(0, as) = 1. (54) 

It follows that the fixed order expansion of x niust be poorly behaved near M = 
0: a simple model of this behaviour is to think of replacing as/M with as/{M+ 
as) = as/M — o?^jyP' + ... in order to satisfy the momentum conservation 
constraint. 

We thus propose a reorganization of the expansion of x i'^to a "double 
leading" (DL) expansion, organized in terms of "envelopes" of the contributions 
summarized in fig. 3b: each order contains a "vertical" sequence of terms of 
fixed order in as , supplemented by a "diagonal" resummation of singular terms 
of the same order in as if as/M is considered fixed. To NLO the new expansion 
is thus 

X(M, as) = [asXoiM) + Xs (ff ) - ^] 



+as 



asXi{M) + xss {^) ~ as + %) - fo] + ■ ■ ■ (55) 



where the LO and NLO terms are contained in the respective square brackets. 
Thus the LO term contains three contributions: xoiM) is the leading BFKL 
function eq. (p^), Xs{as/M) eq. are resummed leading singularities de- 
duced from the one loop anomalous dimension, and ricas/ {'kM) is subtracted to 



avoid double counting. At LO the momentum conservation constraint eq. (54) 
is satisfied exactly because 7o(l) = and [xo(-^^) — ^ M"^ near M — Q 
(see eq.(p^)). At NLO there are again three types of contributions: xi(^) 



from the NLO fixed order calculation (eq. (26)), the resummed subleading sin- 
gularities Xss{aslM) deduced from the two loop anomalous dimension, and 
three double counting terms, /q = fi — — n/(13 -f 10n^)/(367r^n^) and 
/2 = nllW -\- 2n//n|?)/(127r^) (corresponding to those terms with {m,n) = 
(1, 0), (2, 1), (2, 2) respectively in fig. 3b). Note that at the next-to-leading 
level the momentum conservation constraint is not exactly satisfied because 
the constant contribution to xi does not vanish in MS, even though it is nu- 
merically very small (see fig. 4). It could be made exactly zero by a refinement 
of the double counting subtraction but we leave further discussion of this point 
for later. 

Plots of the various LO and NLO approximations to x are shown in fig. 4. 
In this and other plots in this paper we take as = 0.2, which is a typical value 
in the HERA region, and the number of active flavours n/ = 4. We see that, as 
discussed above, the usual fixed order expansion eq. ( p^ in terms of xo and xi 
is very unstable. However, the new expansion eq. ( ^5| ) is stable up to Af ^ 0.3 — 
0.4. Furthermore, in this region, x evaluated in the double leading expansion 
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Figure 4: Plots of different approximations to the v function discussed in the text: the BFKL 
leading and next-to-leading order functions eq. (|l9[), OsXO ^i^d OsXO + <^sXi (dashed); the 
LO and NLO dual cisXs and olsXs + CfgXss of the one and two loop anomalous dimensions 
(solid), and the double-leading functions at LO and NLO defined in eq. ( ^^ (dotdashed). 
Note that the double-leading curves coincide with the resummed ones at small M, and with 
the fixed order ones at large M. 



(pq) is very close to the resummations of leading and subleading singularities 
) obtained by duality eq. (|5l| , |5^ ) from the one and two loop anomalous 
dimensions. This shows that in this region the dominant contribution to and 
thus to 7 , comes from the resummation of logarithms of Q^//!^ with 3> /i^. 

Beyond M ~ 0.4, the size of the contributions from collinear singular 
and nonsingular terms becomes comparable (after all here ^ fi^), but the 
calculation of the latter (from the fixed expansion eq. (|l^)) has become un- 
stable due to the influence of the singularities at M — 1. No complete and 
reliable description of x seems possible without some sort of stabilization of 
these singularities. However, since they correspond to infrared singularities of 
the BFKL kernel (specifically logarithms of Q^/fi^ with <^ fi^) this would 
necessarily be model dependent. In particular, such a stabilization cannot be 
easily deduced from the resummation of the M = singularities: the origi- 
nal symmetry of the gluon-gluon amplitude at large s is spoiled by running 
coupling effects and by unknown effects from the coefficient function through 
which it is related to the deep-inelastic structure functions, in a way which 
is very difficult to control near the photoproduction limit M ^ 1. We thus 
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prefer not to enter into this problem: rather we wiU discuss later a practical 
procedure to bypass it. 

The results summarized in fig. 4 clearly illustrate the superiority of the 
new double leading expansion of x over the fixed order expansion, and al- 
ready indicate that the complete x function could after all lead to only small 
departures from ordinary two loop evolution. 



6 Improving the Splitting Function Expansion 

Having constructed a more satisfactory expansion eq. (^) of the kernel x^ 
we now derive from it an improved form of the anomalous dimension 7 to be 
used in the evolution eq. in order to achieve a more complete description of 
scaling violations valid both at large and small x. In principle, this can be done 
by using the duality relation eq. (|33|), which simply gives the function 7 as the 
inverse of the function x- However, in order to derive an analytic expression 
for ^{N,as) which also allows us to clarify the relation-tp previous attempts 
we start from the naive double-leading expansion of 7 in which terms are 
organized into "envelopes" of the contributions summarized in fig. 3a in an 
analogous way to the double leading expansion ( |55| ) of x'- 

7(iV,a.) = K7o(A^)+7.(t)- 

+as [asii{N) + 7s. (t) - + - cq] + • • • , (56) 



where now 62 = .92*'' = 0, ei = gj'*"-' — nfadb + 13/(2n^))/(187r^) and eo = 
— {^n^+n f) / {Girn'^) . In this equation, the leading and subleading singularities 
7s and 7^5 are obtained using duality eq. (^) from xo ^-ud Xij ^ud summed 
up to give expressions which are exact at NLLx. These are then added to the 
usual one and two loop contributions, and the subtractions take care of the 
double counting of singular terms. 



It can be shown that the dual of the double leading expansion of x eq. (55) 
coincides with this double leading expansion of 7 eq. ( p6| ) order by order in 
perturbation theory, up to terms which are higher order in the sense of the 
double leading expansions. However, it is clear that these additional sublead- 
ing terms must be numerically important. Indeed, it is well know that at 
small N the anomalous dimension in the small- a; expansion eq. ( ^0|) is com- 
pletely dominated by ^ss{(^s/N) which grows very larga,and negative, leading 
to completely unphysical results in the HERA region EZl. It is clear that this 
perturbative instability will also be a problem in the double leading expansion 
eq. (|56|). On the other hand, we know from fig. 4 that the exact dual of x in 
double leading expansion is stable, and not too far from the usual two loop 
result. The origin of this instability problem, and a suitable reorganization 
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of the perturbative expansion which allows the resummatipn of the dominant 
part of the subleading terms have been discussed in ref. EJ. After this resum- 
mation, the resulting expression for 7 in double leading expansion will be very 
close to the exact dual of the corresponding expansion of x- 

To understand this point we consider the asymptotic behaviour of and 
Pss at large ^. Starting from the definition of the splitting function as the 
inverse Mellin transform of the anomalous dimension 7 as given in eq.(|4^), we 
have: 

where we have used eq.(p3|) to change the integration variable according to 
dN = asXo(7s)fi7s- Similarly for Pg^, by using the form of 7ss given in eq.(^), 
we obtain: 

PssiaM^ ^^''''rssijt) = + |5e->^»('^=)«xi(7.) (58) 

The asymptotic behaviour of Ps and Pss at large ^ can now be obtained by 
applying the saddle point method. Remembering that the only real minimum 
of xoils) in the range < 7^ < 1 is at 7s = ^, and expanding around the 
minimum according to 

IsXoils) = ixo(5)(7s-^) + Xoih) (7.-5)'+--; Xihs) = Xiih) + -- (59) 
one finds (see the formulae on the saddle point method in Appendix): 

PsiasO ~ e"=>^«(5)« =^->o. ^ ^ asxi{l)t (60) 

We see that at sufficiently small values of x Pss overwhelms Ps causing the 
perturbative expansion to become unstable. Actually this occurs at not so 
small values of x because xi/xo is so large. j— . 

On the basis of these results the procedure of ref. c3 can be interpreted 
in a simple way whenever the all-order "true" function x{M^as) possesses a 
minimum at a real value of M, Mmin, with < Mmin < 1 (although the final 
result for the anomalous dimension will retain its validity even in the absence 
of such minimum). Using A to denote this minimum value of x, 

A = x(M™„,a,) = Ao + AA, Aq = 0^X0(5) = In 2. (61) 

The instability turns out to be due to the fact that higher order contri- 
butions to 7 must change the asymptotic small x behaviour from x"^" to 
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^.-A ^ ^-\„^A\i ^ + AA^ + ....], with AA = a^xid) H .Having 

understood the source of the problem we can now cure it by a suitable re- 
summation procedure. In terms of splitting functions the proposed resummed 
expansion is simply 

xP{x, as) = ase^^^[Ps{asO + + • ■ •] 

= a,e«^^[P,(a,e) + ^sPssiasO - CAXP^ia^O + •■•]■ (62) 

The expansion is now stable 0, in the sense that it may be shown that 
PssictsO / -Ps{(^sO remains bounded as ^ — > oo: subleading corrections will 
then be small provided only that as is sufficiently small. Equivalently, this 
procedure consists of absorbing the value of the correction to the value of x at 
the minimum into the leading order term in the expansion of x- 

x{M,as) = asXo{M) + alxi{M) + . . . 

= (asXo(M) + AA)+a2xi(M) + .--, (63) 

where Xn{M) = Xn{M) — Cn, with c„ chosen so that Xn(Af) no longer leads 
to an 0(a") shift in the minimum. Since the position Mmin of the all-order 
minimum is not known, one must in practice expand it in powers of as around 
the leading order value M = so at higher orders the expressions for the 
subtraction constants c„ can become quite complicated functions of Xi a-nd 
their derivatives at M — However at NLO we have simply ci — XiC^)? so 
AA = a2xi(i)-l----. 

The corresponding expansion of 7 in resummed leading and subleading 
singularities can also be obtained from the duality eqs. (^ , ^ ,. . . ) by treating 
Xo + AA as the LO contribution to x^ and the subsequent terms Xi as perturba- 
tive corrections to it. Of course, since the reorganization cq. ( |6^ ) amounts to a 
reshuffling of perturbative orders, to any finite order the anomalous dimension 
obtained in this way will be equal to the old one up to formally subleading 
corrections. Explicitly, we find in place of the previous expansion in sums of 
singularities eq. (^) the resummed expansion 

j{N,as)^js{j^)+asjss{jr^)+..., (64) 

where 

(^) - 7. (^) - .■'fl (65) 

The shift in the denominators from N to N — AA results from combining the 
exponential e^^^ of the Mellin transform with the exponential e^'*'^ factored 
out in eq. (p2|) . 
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We can thus replace the unresummcd singularities 7^ and jss in eq. (56) 
with the resummed singularities eq. ( |64[ ) to obtain a double leading expansion 
with stable small x behaviour: 



JV-AA 



eo 



(66) 



Momentum conservation is violated by the resummation because 7s and 7^3 
and the subtraction terms do not vanish at iV = 1. It can be restored by 
simply adding to the constant eo a further series of constant terms beginning 
at 0{al): these are all formally subleading in the double leading expansion. 
This constant shift in 7 is precisely analogous to the shift made on x in eq. (61) 
which generated the resummation. 

It is important to recognize that there is inevitably an ambiguity in the 
double counting subtraction terms in cq. (|6^). For example, at the leading 

we could 



ttN 



order of the double leading expansion instead of subtracting 
have subtracted -;f(jfr^x)' since this differs only by formally subleading terms: 
AA = 0(a2), so 



N 



N - AX 



1 



AA 



TV- AA 



(67) 



Following the same type of subtraction at NLO, the resummed double leading 
anomalous dimension may thus be written as 



7(iV,a,)= a,^o(iV) + 7,(^)-«,^(5^i-I) 



+as 



«.7i(iV)+7..(iv^) + ;r(#^ 



(68) 



where the new contribution proportional to ricAA in the NLO subtraction is 
due to the fact that to NLO accuracy we cannot simply replace asjN with 
as/{N — AA), but we must retain the NLO correction term in eq. (pT]). The 
characteristic feature of this alternative resummation is that the fixed order 
anomalous dimensions 70, 71 are preserved in their entirety, including the posi- 
tion of their singularities. As with the previous expansion eq. (^6|) momentum 
conservation may be imposed by adding to eo a series of terms constant in N 
and starting at 0{a^). 

This completes our procedure of inclusion of the most important part of 
the subleading corrections, as we shall see shortly by a direct comparison of the 
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Figure 5: Comparison of the anomalous dimension 7 evaluated at NLO in the resummed 
expansion eq. (fe™ for three different values of A (dashed) with the usual fixed order per- 
turbative anomalous dimension (also at NLO) eg. (^) (dotted) and that obtained by exact 
duality from x NLO in the expansion eq. (pSl) as displayed in fig. 4 (solid). The unre- 
summed 7 eq. ( p^ ) is also shown at NLO. Notice that the A = 0.21 curve is very close to the 
two loop anomalous dimension down to the branch point at A'^ = A. 

resummed expansions eq. ( p6| ) and eq. (|6^) with the exact dual of x evaluated 
according to eq. (|55|). In the sequel we will discuss the phenomenology based 
on the two resummed expansions eq. (^) and eq. ( |69| ) on an equal footing, 
taking the spread of the results as an indication of the residual ambiguity 
due to subleading terms. Although formally the differences between the two 
expansions are subleading, we will find that in practice they may be quite 
substantial, because AA may be large. 

7 Application to Phenomenology 

So far we have constructed resummations of the anomalous dimension and 
splitting function which satisfy the elementary requirements of perturbative 
stability and momentum conservation. This construction relies necessarily on 
the value A of % near its minimum, since it is this which determines the small 
X behaviour of successive approximations to the splitting function. In order to 
obtain a formulation that can be of practical use for actual phenomenology, we 
will need however to improve the description of x(-^^) in the "central region" 
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near its minimum Mmin, since as we already observed, we cannot reliably 
determine the position and value of the minimum of x without a stabilization of 
the M = 1 singularity. Indeed, we can see from fig. 4 that in the central region 
X evaluated in the double leading expansion is dominated by the presumably 
unphysical M = 1 poles of Xi E^nd at NLO this means that it actually has no 
minimum, becoming rapidly negative. However, one can use the value A of 
the true x the minimum as a useful parameter for an effective description 
of the X function around M = 1/2. Indeed, AA as estimated from its next- 
to-leading order value agXi(l/2) turns out to be of the same order as Aq for 
plausible values of a^, a feature which can be also directly seen from fig. 4. 
This supports the idea that A and AA are not truly perturbative quantities: in 
general we expect that the overall shift of the minimum will still be of the order 
of Ao and negative. It is this order transmutation that makes the impact of the 
resummations eq. ( |6^ , |69| ), and the differences between them, quite substantial. 
After the reinterpretation of A as a parametrization of our ignorance of x in the 
region near AI — i the resulting x function is fixed near M = by momentum 
conservation and near the central region by the value of A. Thus in the region 
of interest for deep inelastic scattering it cannot be much different from the 
true X function. 

In fig. 5 and fig. 6 we display the results for the resummed anomalous 



dimensions in the two different expansions eq. (66) and eq. ( |69| ) respectively, 
each computed at next-to-leading order. In both figures we show for compar- 
ison the fixed order anomalous dimension as7o(iV) -|- a'^"fi{N) eq. (^. Also 
for comparison, we show the exact dual of x computed at NLO in the double 
leading expansion eq. (|55|), obtained from eq. ( |33| ) by exact numerical inver- 
sion. This curve is thus simply the inverse of the corresponding curve already 
shown in fig. 4. 

In fig. 5 we show the anomalous dimension computed at NLO using the 
resummation eq. (|66|), for A = Aq and A — 0. The first value corresponds 
to the LO approximation to A, while the second value is close to the NLO 
approximation when as is in the region as ^ 0.1 — 0.2. We might expect the 
value of A as determined by the actual all-order minimum of x to lie within 
this range. Note that, in general, the resummed anomalous dimension has a 
cut starting at iV = A, which corresponds to the power rise; for this reason 
our plots stop at this value of N. The A = curve, corresponding to the next- 
to-leading order approximation to A, is seen to be very close to the exact dual 
of X at NLO in the expansion eq. (|55|), as already anticipated. This is to be 



contrasted with the corresponding unresummed anomalous dimension eq. (^6|) , 
which is also displayed in fig. 5, and is characterized by the rapid fall at small 
N discussed already in ref.E-OO. This comparison demonstrates that indeed the 
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Figure 6: Comparison of the anomalous dimension 7 evaluated at NLO in the resummed ex- 
pansion eq. ( pg] ) for two different values of A (dashed) with the usual fixed order perturbative 
anomalous dimension (also at NLO) eq. (^) (solid) and that obtained by exact duality from 
X at NLO as in fig. 4 (dotdash). Notice that the A = curve is virtually indistinguishable 
from the fixed order anomalous dimension for all values of N. 



perturbative reorganization eliminates this pathological steep decrease. The 
resummed curve with A = and the exact dual of x become rather different for 
small N ^ 0.2. However, this is precisely the range of N which corresponds 
to the central region of M where we cannot trust the next-to-leading order 
determination of x- Finally, we show that we can choose a value of A ~ 0.21 
such that the resummed anomalous dimension closely reproduces the two loop 
result down to the branch point at = A. This shows that the absence of 
visible deviations from the usual two loop evolution can be accommodated 
by the resummed anomalous dimension. However this is not necessarily the 
best option phenomenologically: perhaps the data could be better fitted by a 
different value of A if a suitable modification of the input parton distributions 
is introduced. It is nevertheless clear that large values of A such as A « Ao 
can be easily excluded within the framework of this resummation, since they 
would lead to sizeable deviations from the standard two loop scaling violations 
in the medium and large x region. 

The splitting functions corresponding to the anomalous dimensions of fig. 5 
are displayed in fig. 7. The basic qualitative features are of course preserved: 
in particular, the curves with small values of A = and A = 0.21 are closest 



24 



to the two loop result. However, on a more quantitative level, it is clear 
that anomalous dimensions which coincide in a certain range of N, but differ 
in other regions (such as very small N) may lead to splitting functions which 
differ over a considerable region in x. In particular, the A = 0.21 curve displays 
the predicted growth at sufficiently large ^ {x ^ 10~^). The dip seen in 
the figure for intermediate values of ^ is necessary in order to compensate 
this growth in such a way that the moments for moderate values of N remain 
unchanged. NoteJdiat the x~^ behaviour of the splitting function at small x is 
corrected by logstiJ: Ps ~ ^~^/^x~'^. If A = this logarithmic drop provides 

the dominant large ^ behaviour which appears in the figure. 

If the anomalous dimensions are instead resummed as in eq. (|69|), the 
results are as shown in fig. 6, again for the two very different values of A, 
A = and A = Aq. When A = the resummed anomalous dimension is now 
essentially indistinguishable from the two loop result. This is due to the fact 
that the simple poles at = which are now retained in 70 and 71 provide 
the dominant small N behaviour. The branch point at iV = A in 7^ and 7<,s 
is then relatively subdominant. This remains of course true for all A < 0, 
and in practice also for small values of A such as A ^ 0.1. When instead 
A = Ao the result does not differ appreciably from the resummed anomalous 
dimension shown in fig. 5, since now the dominant small N behaviour is given 
by the branch point at iV = Ao, which is not affected by changes in the double 
counting prescription. Summarizing, the peculiar feature of the resummation 
eq. (69) is that it leads to results which arc extremely close to usual two loops 
for any value of A < 0, without the need for a fine-tuning of A. 

Finally, in fig. 8 we display the splitting functions obtained from the re- 
summed anomalous dimensions of fig. 6. The A = Aq case is, as expected, 
very close to the corresponding curve in fig. 7. However the A = curve is 
now in significantly better agreement with the two loop result than any of the 
resummed splitting functions of fig. 7, even that computed with the optimized 
value A = 0.21. Moreover, this agreement now holds in the entire range of ^. 
This is due to the fact that the corresponding anomalous dimension is now very 
close to the fixed order one for all > 0, and not only for > A = 0.21. This 
demonstrates explicitly that one cannot exclude the possibility that the known 
small X corrections to splitting functions resum to a result which is essentially 
indistinguishable from the two-loop one. This however is only possible if A 0. 



To summarise, we find that the known success of perturbative evolution, 
and in particular double asymptotic scaling at HERA can be accommodated 
within two distinct possibilities, both of which are compatible with our current 
knowledge of anomalous dimensions at small x, and in particular with the 
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Figure 7: The splitting functions corresponding to the anomalous dimensions of fig. 5. 
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Figure 8: The splitting functions corresponding to the anomalous dimensions of fig. 6. 
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inclusion of corrections derived from the BFKL equation to usual perturbative 
evolution. One possibility, embodied by the resummed anomalous dimension 
eq. (69) with A ^ 0, is that double scaling remains a very good approximation 
to perturbative evolution even if the x — s- limit is taken at finite . The 
other option, corresponding to the resummation eq. (66) with a small value 
of A, is that double scaling is a good approximation in a wide region at small 
X, including the HERA region, but eventually substantial deviations from it 
will show up at sufficiently small x. In the latter case, the best- fit parton 
distributions might be significantly differ from those determined at two loops 
even at the edge of the HERA kinematic region. Both resummations are 
however frilly compatible with a smooth matching to Regge theory in the low 
Q"^ regionEil. 



8 Conclusion 

In conclusion, we have presented a procedure for the systematic improvement of 
the splitting functions at small x which overcomes the difficulties of a straight- 
forward implementation of the BFKL approach. The basic ingredients of our 
approach are the following. First, we achieve a stabilization of the pertur- 
bative expansion of the function x near M = through the resummation of 
all the LO and NLO coUinear singularities derived from the known one- and 
two-loop anomalous dimensions. The resulting x function is regular at M = 0, 
and in fact, to a good accuracy, satisfies the requirement imposed by momen- 
tum conservation via duality. Then, we acknowledge that without a similar 
stabilization of the M = 1 singularity it is not possible to obtain an reliable 
determination of x in the central region M ^ 1/2. However, we do not have an 
equally model-independent prescription to achieve this stabilization at M = 1. 
Nevertheless, the behaviour of x in the central region can be effectively pa- 
rameterized in terms of a single parameter A which fixes the asymptotic small 
X behaviour of the singlet parton distribution. This enables us to arrive at an 
analytic expression for the improved splitting function, which is valid both at 
small and large x and is free of perturbative instabilities. 

This formulation can be directly confronted with the data, which ulti- 
mately will provide a determination of A along with cng and the input parton 
densities. The well known agreement of the small x data with the usual 
evolution equations suggests that the optimal value of A will turn out to be 
small, and possibly even negative for the relevant value of a^. Such a value of A 
is theoretically attractive, because it corresponds to a structure function whose 
leading- twist component does not grow as a power of x in the Regge limit: it 
would thus be compatible with unitarity constraints, and with an extension of 
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the region of applicability of perturbation theory towards this limit. 

Several alternative approaches to deal with the same problem through the 
resummation of various classes of formally subleading contributions have been 
recently presented in the literature. Specific xtroposals are based on making a 
particular choice of the renormalization scaleEa, or on a different identification 
of the large logs which, are resummed by the ^ ptpljition equation (|l^) , either by 
a function of ^ itselfEa, or by a function of Q^tSGZl, or bothEa. The main shoct^ 
coming of these approaches is their model dependence. For instance, in ref. EH 
the value of A is calculated, and x is supposedly determined for all < M < 1. 
This however requires the introduction of a symmetrizaiinn of x, which we con- 
sider to be strongly model dependent: indeed, in ref. EZl it is recognized that 
their value of A only signals the limit of applicability of their computation. 
We contrast this situation with the approach to resummation presented here, 
which makes maximal use of all the available model-independent information, 
with a realistic parameterization of the remaining uncertainties. We expect 
further progress to be possible only on the basis of either genuinely nonpertur- 
bative input, or through a substantial extension of the standard perturbative 
domain. 
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Appendix: The Saddle Point Method 

We want to approximately evaluate the integral 

I{s) = / —exp WiNMN) (69) 

as an expansion in the small parameter 1/s, where f{N) and g{N) are given, 
sufficiently well behaved, real functions. Assume that there is one (or more) 
points where f{N) is stationary: f'{No) — 0. Expanding f{N) and g{N) 
around Nq we have 

f{N) - f{No) + ^f"{No){N - Nof + (70) 

giN) = f{No) + g'(iVo) (N - No) + y"{No) (N - iVo)^ + 
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Performing the gaussian integrals one obtains for each iVo a eontribution: 
exp5/(iVo) . giNo) , .g"(7Vo) , , 



lis) 



which gives the desired expansion in powers of 1/s. 
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